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ABSTRACT 

The multi-phase interstellar medium (ISM) is highly structured, on scales from the 
size of the Solar System to that of a Galaxy. In particular small-scale structures are 
difficult to study and hence are poorly understood. We used the multiplex capabilities 
of the AAOmega spectrograph at the Anglo-Australian Telescope to create a half- 
square-degree map of the neutral and low-ionized ISM in front of the nearby (~ 5 
kpc), most massive Galactic globular cluster, a; Centauri. Its redshifted, metal-poor 
and hot horizontal branch stars probe the medium-strong Can K and Nai D2 line 
absorption, and weak absorption in the A5780 and A5797 Diffuse Interstellar Bands 
(DIBs), on scales around a parsec. The kinematical and thermodynamical picture 
emerging from these data is that we predominantly probe the warm neutral medium 
and weakly-ionized medium of the Galactic Disc-Halo interface, ~ 0.3-1 kpc above the 
mid-plane. A comparison with Spitzer Space Telescope 24- /jm and DIRBE/IRAS maps 
of the warm and cold dust emission confirms that both Na 1 and Ca II trace the overall 
column density of the warm neutral and weakly-ionized medium. Clear signatures are 
seen of the depletion of calcium atoms from the gas phase into dust grains. Curiously, 
the coarse DIRBE/IRAS map is a more reliable representation of the relative reddening 
between sightlines than the Nai and Can absorption- line measurements, most likely 
because the latter are sensitive to fluctuations in the local ionization conditions. The 
behaviour of the DIBs is consistent with the A5780 band being stronger than the 
A5797 band in regions where the ultraviolet radiation level is relatively high, as in 
the Disc-Halo interface. This region corresponds to a cr-type cloud in which Cai and 
small diatomic molecules such as CH and CN are usually absent. In all, our maps and 
simple analytical model calculations show in unprecedented detail that small-scale 
density and/or ionization structures exist in the extra-planar gas of a spiral galaxy. 

Key words: stars: horizontal branch - ISM: molecules - ISM: structure - Galaxy: 
disc - Galaxy: halo - globular clusters: individual: w Cen (NGC 5139) 



1 INTRODUCTION 

"Overhead the sky was full of strings and shreds of vapour, 
flying, vanishing, reappearing, and turning about an axis 
like tumblers, as the wind hounded them through heaven." 
(Stevenson 1879). The interstellar medium (ISM) is indeed 
an extra-ordinarily dynamic place. Cooling and heating pro- 
cesses are continuously at work to transform gas between 
different states, from hot ionized gas with T ~ 10 6 K visible 
in X-rays, to a cold neutral medium with T < 100 K visi- 
ble in Hi and far-IR emission (Heiles & Troland 2005). The 
hotter media tend to fill space, whereas the colder media 



tend to form discrete structures (Cox 2005). The dynamics 
are important for Galactic recycling such as the re-accretion 
of gas from the Galactic Halo as fountains of hot gas from 
supernovae in the Galactic Disc cool and condense. They 
also influence the formation of small-scale structure in the 
form of filaments and small clumps, which must be transient 
features as they are out of pressure equilibrium with their 
surroundings (Stanimirovic & Heiles 2005). Some of these 
structures have AU scales; we know very little about the 3D 
shapes of these tiny-scale structures, although a sheet-like 
geometry seems likely (Heiles & Troland 2005). 
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Absorption-line studies, e.g., against pulsars in radio or 
stars in optical/ultraviolet, can probe the small-scale struc- 
ture; by using the relative proper motion of the background 
source, ID tiny-scale structure can be traced (Heiles 1997; 
Crawford 2003). Mapping small- and tiny-scale structure di- 
rectly is often limited by large separations between back- 
ground sources, and limited angular resolution. Together 
with mapping high-surface brightness nebulae (e.g., Lazio 
et al. 2009), globular clusters offer some of the best possi- 
bilities for mapping experiments. 

Langer, Prosser & Sneden (1990) observed several sight- 
lines towards the globular clusters M 15 and M 92, finding 
evidence for significant column density and velocity vari- 
ations on scales as small as 0.2 pc. Bates et al. (1995) 
presented similar findings for sightlines towards the glob- 
ular cluster M 13. Variations on even smaller scales, down 
to ~ 0.01 pc, were detected towards globular clusters M4 
(Kemp, Bates & Lyons 1993) and the central parts of M 15 
(Meyer & Lauroesch 1999) and M92 (Andrews, Meyer & 
Lauroesch 2001). The latter team points out that fluctua- 
tions in ionization equilibrium and not total column den- 
sity may be responsible for the observed variations. Points, 
Lauroesch & Meyer (2004) observed 150 sightlines towards 
the binary open cluster h & xP erse ii detecting both small- 
scale fluctuations as well as correlated behaviour over larger 
scales, suggesting variations in the physical conditions within 
sheets of gas. 

We are incredibly fortunate to find ourselves in close 
proximity to uo Centauri, the most massive Galactic globu- 
lar cluster (M ~ 2.5 x 10 6 M ; van de Ven et al. 2006). 
Few globular clusters are nearer to the Sun than u Cen 
(d = 5 kpc; McDonald et al. 2009), resulting in a large num- 
ber of suitable, bright probes of angular scales between an 
arcsecond and a degree — (sub)parsec scales in the inter- 
vening ISM. We here use its many horizontal branch (HB) 
stars. These stars are ideal probes: they are metal-poor, 
-1.8 <[Fc/H]< -0.8, and hot, T > 7000 K, and thus rel- 
atively clear of photospheric absorption in the interstellar 
lines of low excitation that probe the neutral and low-ionized 
ISM. The high retrograde velocity of ui Cen (ulsr ~ 229 
km s _1 ) moves the photospheric absorption out of the inter- 
stellar column. The proper motions of the stars are known 
too (~ 5 mas yr _1 , van Leeuwen et al. 2000), so repeat 
observations may be used to probe the tiny-scale structure. 

The modest extinction, E(B - V) ~ 0.08 mag (McDon- 
ald et al. 2009), and 15° Galactic latitude of u) Cen make the 
interstellar path towards u> Cen relatively simple. Yet several 
components have been identified, from the Local Bubble to 
a more distant spiral arm (Wood & Bates 1994). The in- 
teresting possibility is that some of the material we would 
probe could be in the Disc-Halo interface or associated with 
wCen itself; emission is seen at velocities of ~ 195 km s -1 
in Hi (McDonald 2009; McDonald et al, in prep.) and CO 
(Origlia et al. 1997). Absorption from the ISM in the 3968.47 
& 3933.68 A Can H & K doublet was easily detected in 2dF 
spectra (van Loon et al. 2007)Q. The same would be expected 
for the 5895.92 & 5889.95 Nal Di & D 2 doublet. The 2dF 
spectra were of too low a resolution to resolve the interstel- 



1 Jones (1968) already noted the appearance of interstellar K-linc 
absorption in the spectrum of the F-type giant uj Cen VI. 



lar lines. We have since used its successor, AAOmega, at 
R ~ 8000, to measure ISM features in 452 hot HB stars in 
u Cen, of which we present the results here. 



2 OBSERVATIONS 

2.1 Spectroscopy with AAOmega 

The observations were undertaken, in service mode, at 
the Anglo- Australian Telescope (AAT) on the night of 28 
March 2008. The AAOmega spectrograph was used with the 
multiple-object-spectroscopy fibre feed from the 2dF fibre 
positioner (Saunders et al. 2004; Sharp et al. 2006). 

The 5700 A AAOmega dichroic was used. The red and 
blue arms of AAOmega were used with the 3200B and 2000R 
VPH gratings centered and blazed at 4040 A and 5990 A 
yielding spectra in the 3890-4170 A and 5740-6220 A re- 
gions at spectral resolutions (per 3.4-pixel resolution ele- 
ment) of R ~ 8200 (Av ~ 37 km s _1 ) and R ~ 7245 
(Ai> ~ 41 km s~ ), respectively. The exact ranges of individ- 
ual spectra differ. The Ca n H & K lines, which were the pri- 
mary targets for the blue arm observations, were placed off 
centre on the CCD to avoid regions of poor cosmetic perfor- 
mance (bad columns). Skies were clear with seeing varying 
during the night from 1.4" to 2.0"; the 2dF fibre entrances 
subtend 2.0" on the sky. 

Two fibre configurations were employed. Each observing 
block consists of a flat field frame (quartz-halogen lamp), an 
arc frame for wavelength calibration (CuAr, FeAr, He and 
Ne), a set of twilight sky flat field frames (to normalize the 
relative fibre transmissions for sky subtraction) and a series 
of 1800 second science frames. The total observing times for 
the two configurations sampling different stars were 6 x 1800 
seconds and 4 x 1800 + 1 x 1200 seconds, respectively. 

The data were processed using the 2dfdr data reduc- 
tion package. Two additional steps were performed under 
IRAF prior to processing with 2dfdr. Firstly, a 2D bias frame 
was subtracted to correct a number of poor but recoverable 
columns in the data set. The 2D bias frame was created for 
this purpose from a stack of 30 bias frames via a cr-clipped 
mean rejection. Secondly, the iraf fixpix task was used to 
interpolate across a single bad column at the wavelength 
of the Can K line. This isolated bad column affected the 
Call K line in around 25% of spectra, with spectral curva- 
ture across the CCD moving it clear for the remaining stars. 
An updated bad-pixel mask was also created for 2dfdr to 
account for the effects of these modifications. 

2dfdr performs the standard reduction operations: 
overscan correction; fibre trace and extraction; flat fielding 
and wavelength calibration. Cosmic ray rejection was ap- 
plied to the 2D science frames by 2dfdr prior to extraction, 
following the prescription of van Dokkum (2001). In each 
configuration, 25 fibres were devoted to blank sky regions to 
create an average spectrum of the sky. Sky subtraction was 
performed using this sky spectrum. The relative intensity 
of the sky to subtract was determined from twilight flat- 
field frames obtained at the end of the night. Science frames 
were combined using a relative flux weighting derived from 
spectral intensities in each frame. 

The spectra were not Doppler-corrected, i.e. they reflect 
the observer's motion through space. All velocities used in 
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Table 1. Catalogue of observed stars in ui Cen. Only the first eight rows are displayed; the full catalogue of 452 objects is made available 
at CDS. The first columns list the stellar probe's LEID number (van Leeuwen ct al. 2000) and its equatorial and Galactic coordinates, 
and the reddening E(B-V) from the DIRBE/IRAS maps of Schlegel et al. (1998). Subsequent columns list the central velocity with 
respect to the Local Standard of Rest, line width (ct) and equivalent width (W) of the Can K and Nal D2 interstellar lines, and the 
equivalent widths of the A5780 and A5797 Diffuse Interstellar Bands. 



Star 


Equatorial coordinates 


Galactic coordinates E(B-V) 




Cai K 






Nai D 2 




A5780 


A5797 




a (J2000) 


5 (J2000) 


I 


b 




V 


ct 


W 


V 


CT 


W 


W 


W 


LEID 


(deg) 


(deg) 


(deg) 


(deg) 


(mag) 


(km/s) 


(mA) 


(mA) 


(km/s) 


(mA) 


(mA) 


(mA) 


(mA) 


40004 


201.06673 


-47.45981 


308.669 


15.0463 


0.1301 


-19.4 


28.6 


525 ± 16 


-13.6 


22.7 


474 ± 15 


82 ± 15 




42009 


201.11112 


-47.47200 


308.698 


15.0302 


0.1304 


-32.1 


30.3 


520 ± 18 


-13.2 


23.0 


418 ± 12 


59 ± 24 




31006 


201.12876 


-47.39312 


308.721 


15.1068 


0.1256 


-24.3 


29.6 


521 ± 15 


-13.4 


23.6 


375 ± 12 


75 ± 16 




64010 


201.18389 


-47.63934 


308.726 


14.8578 


0.1417 


-27.9 


32.2 


627 ± 54 


-16.5 


26.2 


538 ± 50 


52 ±22 


37 ± 19 


20006 


201.14240 


-47.30854 


308.743 


15.1894 


0.1227 


-17.8 


27.2 


494 ± 16 


-12.1 


22.6 


393 ± 13 


95 ± 25 


43 ± 15 


38008 


201.18022 


-47.44549 


308.750 


15.0503 


0.1286 


-25.4 


29.9 


586 ± 18 


-17.3 


27.0 


539 ± 11 






26009 


201.19651 


-47.35413 


308.774 


15.1393 


0.1227 


-23.9 


27.7 


517 ± 15 


-14.0 


25.3 


351 ± 11 


57 ± 17 


25 ± 10 


35011 


201.21149 


-47.42335 


308.775 


15.0694 


0.1272 


-20.6 


29.7 


575 ± 15 


-18.3 


25.4 


559 ± 17 


82 ± 19 





this work and in the catalogue have been corrected to the 
Local Standard of Rest (LSR), though, by subtracting the 
— 10 km s _1 motion of the observer with respect to the LSR 
at the time of observation. 

The continuum signal-to-noise ratio as measured on the 
extracted spectra is typically 30-50 per pixel around the 
Ca 11 K line and 60-100 per pixel around the Na I D2 line. 

2.2 Spectroscopic targets 

Targets were selected from the proper motion catalogue of 
van Leeuwen et al. (2000). The criteria were designed to 
select blue stars, with 14.6 < B < 16.6 mag and —0.2 < 
(B—V) < 0.25 mag, that are likely members of co Cen, with a 
membership probability p > 90% (all stars actually observed 
had p > 95%, and mostly 99 or 100%). A small number of 
targets were rejected due to close neighbours (< 3"), which 
would lead to contaminating light within the 2"-diameter 
fibres. Hence 995 potential targets were available. 

To maximise the uniformity and density of the area on 
the sky covered by our observations, targets were drawn 
from a field with radius R = 0.5° and mutual separation 
> 1'. The bluest stars, with -0.20 < (B — V) < 0.10 mag, 
were given priority, then stars with 0.10 < (B — V) < 0.15 
mag at R > 0.15°, followed by stars with 0.15 < (B - V) < 
0.20 mag at R > 0.20° and stars with 0.20 < (B — V) < 0.25 
mag at R > 0.25°. We used the configure software with 
the simulated annealing method to configure the fields, and 
hence produced a second field with the leftovers from the 
first configuration. In all, 452 targets were observed in two 
field configurations. 

To position the fibres accurately on sky, 6 and 7 fiducial 
stars were observed for the first and second configuration, re- 
spectively. These were chosen from amongst the brighter red 
giants, with 12.8 < V < 13.3 mag and 0.9 < (B — V) < 1.3 
mag, with the same membership criterion as before to guar- 
antee a good relative astrometric and fibre-positioning ac- 
curacy. The sky positions were chosen using the same target 
selection software, avoiding contamination within 3" equiv- 
alent to V = 21 mag in the fibre. 

The observed stars are listed in Table 1, the full version 
of which is made available electronically; their locations in 
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Figure 1. Optical B, B-V colour-magnitude diagram of poten- 
tial targets (small dots) with > 50% membership probability, and 
observed stars (squares). The vertical dotted lines indicate prior- 
ities, with bluer stars given higher priority. 



the B, B-V colour-magnitude diagram and on the sky are 
shown in Figs. 1 and 2, respectively. 



3 RESULTS 

The average spectra in the blue and red regions are shown in 
Fig. 3. Besides the strong Balmer hydrogen absorption lines 
arising in the photospheres of the stars, the interstellar Can 
H & K and Nai Di & D2 lines are very obvious features 
observed with ease in every single spectrum. Two Diffuse 
Interstellar Bands (DIBs) are discernible in the red portion 
of the average spectrum, at 5780 and 5797 A. 
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Figure 2. Positional diagram of potential targets (small dots) 
with > 50% membership probability, 14.6 < B < 16.6 mag and 
—0.2 < (B — V) < 0.25 mag, and observed stars (squares). 
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Figure 3. Overview of the average spectrum, in the blue (top) 
and red (bottom) regions, uncorrected for the spectral through- 
put and response of the atmosphere, telescope and instrument. 
Note that the Can H absorption features overlap the strong pho- 
tospheric He line. 
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Figure 4. Top: The average spectrum around the Call K line. 
The interstellar feature is well separated from the redshifted pho- 
tospheric absorption. The average Gaussian fit is overlaid, demon- 
strating the excellence of the fits. Bottom: The average spectrum 
around the Nal Di and D2 lines. The interstellar D2 line is 
well isolated, but the Di line is affected by the (generally much 
weaker) photospheric D2 absorption. The interstellar line profiles 
are very well matched by single Gaussians. 



3.1 Interstellar Call and Nal absorption 

Average profiles of the Can H & K and Nal Di & D2 lines 
are shown in Fig. 4: the shorter-wavelength interstellar ab- 
sorption lines of these multiplets, Can K and Nai D 2 , are 
well isolated from their photospheric counterparts and are 
seen against a clean photospheric continuum; the weaker 
Call H and Nal Di lines are affected by strong photospheric 
Balmcr He and photospheric Nal D2 line absorption, respec- 
tively. The kinematic offset between stars and ISM is small 
enough to avoid confusion with the 5875.6 A He I D 3 line. 

The equivalent widths, line widths and central velocities 
of the interstellar Ca 11 K and Na I D2 lines were measured by 
fitting a Gaussian (Table 1). This worked satisfactorily, the 
average result of which can be judged in Fig. 4. Note that the 
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Figure 5. Maps of the equivalent width (top row) and mean velocity (bottom row) of the interstellar Call K (left) and Nal D2 (middle) 
lines, and their ratio or difference, respectively (right). The levels vary between ±2 standard deviations around the median (high values 
of absorption equivalent width and hence column densities are light-coloured). The positions of the stellar probes are plotted too. 



line width is given in terms of the a value of the Gaussian, 
and hence is related to the Full- Width at Half Maximum as 
FWHM = 2(2 In 2) 1/2 cr = 2.355 a. This includes the contri- 
bution of the instrumental broadening, cri ~ 15.5 and ~ 17.6 
km s _1 for the Call K and Nal D2 lines, respectively. The 
formal errors on the velocities are a tenth of the spectral 
resolution element, or ~ 4 km s _1 (0). The errors on the 
equivalent widths were derived from the standard deviation 
of the residuals of the Gaussian fitting, summed in quadra- 
ture and weighted by the Gaussian fit. Typical errors are 
-20-30 mA for the Call K line and ~ 15-20 mA for the Nal 
D2 line. 

The equivalent widths in ~ 450 sightlines were then 
used to construct maps of the interstellar absorption as well 



2 A check on the photospheric components reveals coincidence 
between Call and Nal lines in uCen stars within a few km s — 1 . 



as the line ratio (Fig. 5). The equivalent widths of both the 
Can K and Nal D2 lines varied very significantly compared 
to the precision of the measurements, but by less than a 
factor two across the maps (Fig. 5a and b): (WcaiiK) = 
536 mA (range 383-709 mA); (W N alD2) = 500 mA (range 
350-682 mA). The standard deviations in the distributions 
of VKcaiiK and t¥NaiD2 values are 48 and 63 mA, or 9% 
and 13%, respectively, compared to the mean errors in the 
equivalent width values of 25 and 15 mA, or 5% and 3%, 
respectively. The real fluctuations in the Call K and Nal 
D2 maps thus amount to ~ 7% and 12%, respectively. 

The equivalent width values can be converted to column 
densities in the lower energy level of the transition, following 
van Dishoeck & Black (1989): 

17 W 

iV = 1.13 x 10 17 x — , (1) 
where the wavelength is in A, the equivalent width is in 
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mA, and the oscillator strengths are / = 0.65 and 0.67 for 
the Call K and Nal D2 lines, respectively. We thus obtain 
(iVcaiiK) w 6.0 x 10 12 cm" 2 ; <iV N aiD 2 ) « 2.4 x 10 12 cm" 2 . 
The conversion to column density assumes that the lines are 
weak enough to be on the linear part of the curve-of-growth. 
Given the moderate values for the equivalent width and red- 
dening, E(B-V), this is not a bad assumption (see Crawford 
1992a for examples of resolved line profiles at different E(B- 
V) values). It also assumes a single kinematic component, 
whereas in our case the absorption is a blend of at least two 
kinematically distinct components (see below). 

Structure is visible on all scales, between neighbour- 
ing sightlines (but not always) as well as across the 
0.7° field. The atomic line ratio varied slightly more 
(Fig. 5c), <VKcaiiK/^NaiD2> = 1.08 (range 0.72-1.54), or 
(-/VcaiiK/ / j'vNaiD2) ~ 2.5 (range 1.7-3.6). Indeed, although 
there are regions where the two lines vary in concert (e.g., 
the lower- left area and upper border in Fig. 5), there are 
regions where either the Can K (e.g., the top-left corner 
and right border in Fig. 5) or the Nal D2 (Fig. 5: around 
I = 309.25°, b = 15.0°) absorption is relatively strong. 

The Nal D2 line is broader in wavelength than the Ca n 
K line, but not in velocity: (o"NaiD 2 ) = 23 km s" 1 (range 19- 
30 km s" 1 ), whereas (crcaiiK ) = 30 km s" 1 (range 24-37 
km s" 1 ). Interestingly, the ratio of equivalent widths corre- 
lates significantly with that of the line widths (Fig. 6a, with 
Pearson's values r — 0.37, p < 0.0001), and a corresponding 
shift in the central velocity between the two lines is seen 
(Fig. 6b, r = -0.31, p < 0.0001). This suggests that at 
least some sightlines pass through multiple absorption com- 
ponents (clouds), with Na° clouds appearing at somewhat 
more positive velocities, and Ca + clouds at more negative 
velocities. We will come back to this observation in the dis- 
cussion section. 

3.2 Diffuse Interstellar Bands 

The first DIB features to have been found (Heger 1922) 
and shown to be of interstellar origin (Merrill 1934) are at 
5780 and 5797 A. The modest reddening towards our stellar 
probes not surprisingly results in rather weak DIB features, 
but in the average spectrum they are clearly detected at 4 
and 2% depth with respect to the continuum (Fig. 7). 

Gaussian fitting was on the whole still adequate, with 
typical errors of ~20-40 mA for the A5780 DIB and ~15- 
25 mA for the weaker but sharper A5797 DIB. Only fits of 
reasonable fidelity were kept by imposing a threshold on the 
equivalent width, > 20 mA, and constraining the bandwidth 
(a) to within 0.15-1.8 A for the A5780 band and 0.15-1 A for 
the A5797 band. Thus we kept the majority (415) of A5780 
fits but only a third (164) of the A5797 fits. For these, the 
equivalent widths were (VKdibstso) = 100 mA, up to a max- 
imum of 251 mA, with a bandwidth (ctdib578o) = 0.95 A, 
and (Worn 5797) = 33 mA, up to a maximum of 128 mA, 
with a bandwidth (<jdib5797 ) = 0.49 A. The shape of the 
A5797 band deviates from that of a Gaussian profile, not 
least due to a broader, shallower band at 5795 A. Consid- 
ering the weakness of the A5797 band in our spectra the 
Gaussian fitting is as accurate as any method. 

There are substantial variations of the DIB strength 
across the 0.7° field and these are mapped in the stronger 
A5780 DIB in Fig. 8. This DIB is generally stronger where 
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Figure 6. Comparison between interstellar Can K and Nal D2 
lines: equivalent widths (a) and velocities (b) versus the a of the 
Gaussian fits expressed in velocity. 
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Figure 7. The average spectrum around the A5780 and A5797 
DIBs. The average Gaussian fit is overlain. 
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Figure 8. Map of the equivalent width of the A5780 DIB. The 
intensity levels vary between ±2 standard deviations around the 
median. Positions of the stellar probes are plotted only for sight- 
lines with log(WDIB5797/WblB578n) > -0.2. 



the Call K line is weaker (compare Fig. 5). However, some 
strong A5780 absorption is seen in the upper left where the 
Nal D2 is weak but the Call K line is not weak. Despite 
fewer sightlines the DIB map appears to have more small- 
scale structure, possibly because DIB absorption is gener- 
ally associated with the neutral cloud medium which has a 
smaller filling factor than the warmer inter-cloud medium. 
However, much of the scatter is due to the uncertainties 
in the equivalent width measurements. The mean error in 
Worn 5780 is 32 mA, whilst the standard deviation of the 
distribution of Wdib5780 values (i.e. the fluctuations in the 
map) is 35 mA. Given the mean value for Wdib5780 of 100 
mA, this would suggest that the real fluctuations amount to 
"only" ~ 15 mA(15%), still comparable to — or exceeding 
— the fluctuations in Call K and Nal D2. 

A significant positive correlation exists for the ratio 
of A5780 DIB and Nal D2 equivalent widths and the ra- 
tio of Call K and Nal D2 equivalent widths (Fig. 9b, 
r = +0.26, p < 0.0001; very small adjustments were made to 
the r- values to account for correlated errors in WcaiiK and 
W / NaiD2). This suggests that this DIB is strong where Ca + is 
more abundant. There is no significant correlation between 
the ratio of A5780 DIB and Ca n K equivalent widths and the 
Can K/Nal D 2 ratio (Fig. 9a, r = -0.06, p = 0.24). There 
is no clear trend of the DIB ratio with Call K/Nal D2 ra- 
tio either (Fig. 9c) , but if the strongest £-type absorption is 
selected (see below), log(WbiB5797/WDiB578u) > —0.2, then 
these sightlines tend to be where the A5780 DIB is weak, 
but both Call K and Nal D2 absorption are fairly strong 
(Figs. 5 & 8). 
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Figure 9. Comparison of the equivalent widths of the A5780 and 
A5797 DIBs, and the interstellar Call K and Nal D2 lines. 



The Wdib 5797 /Wdib 5780 ratio is an indicator of the ion- 
ization conditions as inferred from the atomic to total (H + 
H 2 ) hydrogen ratio (Cami et al. 1997). The A5780 DIB is 
relatively weak in £-Per type clouds and strong in cr-Sco 
type clouds (Krelowski & Sneden 1995), so the Disc-Halo 
interface region appears to be of a type. The charge state 
of the A5780 and A5797 band carriers is not known, but 
if both are, for example, singly charged cations, then the 
low-UV (() c.f. high-UV (a) classification is consistent with 
the parent neutrals having low (A5797) and relatively high 
(A5780) first photo-ionization energies such that the A5780 
absorption appears when the UV field is strong. The <r-type 
character observed towards ui Cen is consistent with other 
observations of the A5797 and A5780 diffuse bands along 
sightlines out of the Galactic plane. Examples of observa- 
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tions of Milky Way diffuse band absorption are those seen 
towards the Magellanic Clouds (Ehrenfreund et al. 2002; 
Welty et al. 2006), and also SN 1987A (LMC) (Vladilo et al. 
1987) and SN 1986G (NGC5128) (D'Odorico et al. 1989) for 
which the reported Wdib 5797/ Wdib 5780 ratios are 1/4 and 
1/2, respectively. 

The A6196 and A6203 DIBs are also seen in some of the 
spectra. However, this is close to the edge of our spectral 
range, and is in fact not covered in all spectra. In the spectra 
that do cover these bands, they are generally too faint to 
allow a meaningful analysis. 



3.3 Evidence for gas associated with u Centauri? 

The spectral resolution element exceeds the internal velocity 
dispersion of the ui Cen cluster (e.g., van Loon et al. 2007), 
and the generally weak but present photospheric absorption 
in both the Call H and K lines and the Nal Di and D2 lines 
further limits searches for absorption components arising in 
the intra-cluster medium (ICM) within u> Cen. One could 
only hope to find very strong absorption in particularly hot 
stars, or, more sensitively, significantly displaced absorption 
because of the stellar peculiar velocity with respect to the 
cluster's systemic velocity. 

The average spectrum does not reveal any absorption 
in addition to the local interstellar components around the 
Sun's velocity and photospheric components of the stellar 
probes in wCen (Fig. 4). If large variations were to occur 
in the ICM, then there might be a possibility to discover 
ICM components in individual spectra. However, the signal- 
to-noise ratio of the individual spectra is limited, and the 
resulting column density sensitivity is N ij 10 11 cm~ 2 for the 
Can H and K and Nai Di and D2 lines. No evidence was 
found for ICM at a level exceeding this. Assuming a uniform 
column density across a circular area of R — 0.3° radius, 
this corresponds to a limit of Micm Ss 4 M unless Na° and 
Ca + were to be grossly under-representative species. This 
estimate is consistent with the Micm 3 M from Spitzer 
maps of w Cen (Boyer et al. 2008; McDonald et al. 2009), 
and reinforces evidence for efficient clearing of the cluster 
by interaction with the hot Halo (van Loon et al. 2009) . 



4 DISCUSSION 

4.1 Density, ionization, or dust depletion? 

Are we probing density variations, or fluctuations in ion- 
ization and excitation equilibrium, e.g., due to variations 
in the interstellar radiation field? The Solar abundances of 
sodium and calcium are nearly identical, and hence the rel- 
ative column densities of Ca + and Na° readily indicate the 
dust depletion and excitation/ionization conditions. Purely 
thermally-broadened lines of sodium and calcium are unre- 
solved in our spectra. With a mean atomic weight of 23 and 
40, respectively, the line width would be a < 19 km s _1 
even for a 10 6 K plasma. Resolution of these lines therefore 
signifies the presence of multiple clouds along the sightline, 
rather than providing a means to estimate the gas kinetic 
temperature. Furthermore, the ionization state is likely to 
be dominated by the local strength of the interstellar radi- 



ation field, and thus electron density, neither of which are 
particularly well-constrained. 

One may assume that the occupancies of the 3s and 
4s states for the Nai D and Can doublets are mostly re- 
lated to the ionization state rather than a different degree 
of excitation of a particular atomic species. The relevant ion- 
ization potentials are as follows: Ca°, 6.11 eV; Ca+, 11.87 
eV; Ca++, 50.91 eV; Na°, 5.14 eV; Na+, 47.29 eV. This 
means that in the warm ionized medium (significant energy 
density above 13.6 eV), neither Na° nor Ca + are the domi- 
nant species. Still, these lower ionization stages remain well 
populated and, given the high ionization potentials of the 
next higher stage (in both cases), should be fairly stable in 
their occupancy. In the cold neutral medium (with H2), Na 
is dominant but there is no Ca + . Therefore, no correspon- 
dence between Ca 11 and Na 1 can be found in the cold neu- 
tral medium. In the warm neutral/weakly-ionized medium 
(atomic H), Ca + is abundant, but perhaps more suscepti- 
ble to variations in radiation field than Na° given that both 
ionizations to Ca ++ and recombinations to Ca° are frequent 
(c.f. Crawford 1992b). 

The strong Nai 8183+8195 A doublet corresponds to 
the transition upwards from the next state, 3p, but our spec- 
tra do not include this line. The Ca I 4227 A line corresponds 
to the equivalent transition to the Call doublet but in the 
neutral phase, but this line is not covered by our spectra ei- 
ther. The 4227 A line is extremely weak or absent in earlier 
2dF spectra (Figs. 6 & 11 in van Loon 2002). Careful re- 
examination of the 2dF spectra of the few hundred hottest 
stars in van Loon et al. (2007) resulted in a detection rate 
of a few per cent at most, with equivalent widths not ex- 
ceeding a few per cent of that of the Ca 11 K line. Thus we 
have no direct means of constraining the Boltzmann-Saha 
distributions over the electronic energy levels. 

4-1.1 Sodium and calcium as tracers of dust 

Calcium is a highly refractive element and thus easily de- 
pleted from the gas phase onto dust grains, even in diffuse 
clouds. Sodium is much less affected, which allows a test to 
be made by comparing the Na I and Ca 11 equivalent widths 
with the reddening, E(B-V). 

First wc compare our absorption-line maps with the 
partly-overlapping Spitzer Space Telescope 24-fim map (Fig. 
10; c.f. Boyer et al. 2008). The latter shows starlight from the 
cluster in the central portion of the image, but mostly indi- 
cates where relatively warm dust is located. Focussing on the 
band of low 24-/im surface brightness running from about 
I = 308.9°, b = 15.1° to about I = 309.3°, b = 15.2°, this is 
where the Nai absorption is weak but the WcaiiK/WNaiD2 
ratio is fairly high. This suggests that we see warm material 
with little sign of calcium depletion, as expected in regions 
with little dust. Just below this band, around I = 309.2- 
309.3°, b = 15.0°, lies of patch of dust emission; this is where 
Nai is strong and the Wcaii k/Wn»i D2 ratio reaches the 
lowest values in our maps, suggesting high column densities 
of the neutral medium with strong calcium depletion. 

To better quantify the effects of depletion and test 
how well Nai traces the neutral medium, we correlate the 
equivalent widths with reddening values taken from the 
DIRBE/IRAS maps of cold dust (Schlegel, Finkbeiner & 
Davis 1998). These vary by ~ 20% across the maps, between 
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Figure 10. (a:) Spitzer Space Telescope 24-fim image (c.f. Boyer et al. 2008) which, besides starlight from wCen, shows warm dust. The 
sightlines of the absorption-line experiment are overplotted. (6:) Map of the reddening, derived from DIRBE/IRAS data (Schlegel et al. 
1998). Values range from E(B — V) = 0.12 (white) to 0.144 mag (dark). In both panels the centre of tu Cen is indicated by a cross. 



E(B — V) ~ 0.120-0.144 mag; the reddening values are pre- 
cise to 16% (Schlegel et al. 1998), although the statistical 
spread must be considerably smaller within the confines of 
this small portion on the sky (c.f. Fig. 10b). In spite of the 
vastly different angular resolution, the reddening map (Fig. 
10b) reveals clear parallels with the Spitzer 24-/im map (Fig. 
10a). For instance, there is a dense cloud stretching across 
the lower-right corner, and a general lack of dust in the top 
part of the region. Curiously, the reddening map suggests a 
dust cloud very near (but just below) the centre of to Cen. 
This might have been affected by far-IR starlight from the 
cluster, but the Nal absorption (in particular) is also rather 
strong in that patch of sky (Fig. 5). 

The clearest correlation with E(B-V) is for Nal (Fig. 
11a, r — +0.43, p < 0.0001), which increases in strength 
by at least ~ 40%. This confirms that Nal indeed traces 
the column density of the neutral dusty medium (Hobbs 
1974). The fact that it seems to increase more strongly than 
the reddening does, can most simply (and hence plausibly) 
be explained by a need to re-adjust the reddening values 
from the Schlegel et al. (1998) maps towards lower values. 
Indeed, the mean (E{B — V)) ~ 0.13 mag is ~ 0.05 mag 
higher than deduced from accurate analysis of the colour- 
magnitude diagram of stars within uj Cen (McDonald et al. 
2009). If the correction is a simple offset, leaving the spread 
intact, then the reddening would vary by more than 20%, 
possibly by about 40% just as the Nal line strength does. 
Although McDonald et al. (2009) constrain differential red- 
dening gradients across the cluster to AE(B~ V) Si 0.02 mag 
near the cluster centre, Calamida et al. (2005) inferred red- 
dening variations up to a factor almost two from the colours 
of hot HB stars in u> Cen, indicating a clumpy medium; given 
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Figure 11. Comparison of the equivalent widths of the interstel- 
lar Can K and Nal D2 lines as a function of E(B-V) derived from 
DIRBE/IRAS data, with a linear fit to WNalD2 versus E(B-V). 



that the Schlegel et al. maps have a rather coarse resolution 
some of the spread in the Nal line strength may be due to 
such clumps. 

The WcaiiK/WNaiD2 ratio shows a weak anti- 
correlation with E(B-V), at the ~ 20% level (Fig. lib, 
r = —0.24, p < 0.0001), consistent with the effects of dust 
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Figure 12. B, B— V diagram of the spectroscopic targets in 
wCcntauri (top left), and dereddened using the Schlegel et al. 
(1998) DIRBE/IRAS map (top right), or a parameterisation of 
E(B-V) with Nal D 2 (bottom left) or with Call K (bottom right). 



depletion. However, given that the Nal line strengthens even 
more, it means that the Ca n still appears generally stronger 
along dustier sightlines and therefore, to some extent, also 
traces the column density of the neutral medium (c.f. Hobbs 
1974). 

To test whether the Nal or Can lines could be used 
to trace the small-scale distribution of interstellar dust, 
we deredden the B, B-V diagram using the DIRBE/IRAS 
maps, and again using instead a scaling of the reddening 
with the Nal and Can equivalent width (Fig. 12). We pa- 
rameterised E(B - V) = W Na iD2[A]/3.8 and E(B -V) = 
W C aiiK[A]/3.6, and assumed A B = 4.07E(B - V). 

Remarkably, the coarser DIRBE/IRAS map is much 
better in dereddening the colour-magnitude diagram than 
the absorption lines, despite the latter being observed along 
exactly the same sightlines as the stars that are being cor- 
rected. The clearest evidence for this statement is that the 
scatter in the HB sequence increases significantly when using 
the absorption lines, whereas it remains as tight as the origi- 
nal, observed sequence when using the DIRBE/IRAS maps. 
The errors in the equivalent width are typically <, 5%, trans- 
lating into errors in E(B-V) and A B of <, 0.01 and <, 0.05 
mag, respectively, i.e. smaller than the increased scatter. 
The increased scatter, especially when using the Nal absorp- 
tion, would tend to imply that the variations arise from fluc- 
tuations in ionization equilibrium rather than total atomic 
column density, assuming that the atomic column densities 
are well correlated with dust. In any case it confirms that the 
fluctuations in the Can K and Nal D2 lines are a physical 
reality. 
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Figure 13. Comparison between two average spectra, differing 
mainly in the strength of the A5797 DIB. 



4-1.2 The Diffuse Interstellar Bands 

The A5780 and A5797 DIBs are encountered in diffuse 
clouds, with the A5780 DIB appearing stronger where the 
UV radiation field is stronger (cr-type clouds) — which 
presumably leads to the destruction of the carrier of the 
A5797 DIB which is stronger in denser neutral media (£- 
type clouds) (Krelowski & Westerlund 1988; Krelowski et 
al. 1999). If both are singly-charged macromolecules, then 
their different UV threshold for appearance might be related 
to differences in the photo-ionization energy of their parent 
neutral molecules, that of the 5780 A DIB being higher (c.f. 
Cami et al. 1997). Alternatively, it might simply reflect a 
difference in charge, with the A5780 DIB due to a weakly 
ionized or protonated macromolecule and the A5797 DIB 
due to a neutral macromolecule, much like the variations 
suggested in the mid-IR spectra of Polycyclic Aromatic Hy- 
drocarbons (Allamandola, Tielens & Barker 1989). 

The sightlines we probe here would be more translu- 
cent than in the Galactic plane, and hence the a be- 
haviour would be expected to dominate. Indeed, few sight- 
lines in our data have a strong A5797 DIB, consistent with 
a low filling factor of the cold neutral medium. The high- 
est Wdib 5797/WDIB578O ratios are seen where the A5780 
band is weak not because the A5797 band is strong (Fig. 
8). But more typically, it is the A5797 band which varies: in 
Fig. 13 two spectra are compared, showing weak and strong 
A5797 bands for very similar A5780 bands and Ca 11 and Na I 
lines. The former is an average of the spectra of LEID 33137, 
37193, 40185, 41087, 42153, 64042 and 75026, whilst the lat- 
ter is an average of the spectra of LEID 20006, 21022, 26096, 
55048, 60020 and 68026. 

The A5780 DIB shows no correlation with E(B-V) over 
the small range in E(B-V) probed here. Comparing with 
the 24-/xm map (Fig. 10a) there is some correspondence 
between the absence of warm dust around I ~ 309.2°, 
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Figure 14. Correlation between the equivalent width of the 
A5780 DIB and the local maximum gradient in the DIRBE/IRAS 
reddening map. 



b ~ 15.15° and stronger A5780 DIB absorption, and be- 
tween the warm dust emission around I ~ 309.2°, b ~ 14.7- 
14.8° and weaker A5780 DIB absorption. In the latter re- 
gion the A5797 DIB is often seen, confirming the pres- 
ence of a dense medium there. There are about ten times 
more sightlines with Wdib 5797/Wdib5780 < 0.7 than with 
WDIB5797/WDIB5780 > 0.7, which suggests that the denser 
medium traced by the A5797 DIB has a lower filling factor 
than the warmer neutral medium by a factor ~ 10. 

The gradient in reddening can help pinpoint the edges 
of clouds, under the assumption that the reddening reaches 
a maximum in the cloud core and a minimum in the inter- 
cloud regions. The (absolute value of the) gradient would 
be maximum near the boundary between cloud and inter- 
cloud region. There is a very tentative indication that at the 
steepest gradients, dE(B — V r )/dL<;0.06 mag deg -1 , the 
A5780 DIB is stronger than at shallower gradients (Fig. 14): 

(W D IB5780)steep = 111.4±4.6 HlA VerSUS (W D IB5780) shallow = 

96.7 ± 1.8 mA, a difference of 15 ± 5 mA or 3 a. A positive 
correlation, if significant, would be consistent with the A5780 
DIB originating in the warmer envelopes of diffuse clouds. 



4.2 The interstellar path towards u> Centauri 

Much of the structure might be relatively near to us, 
seen against a smoother, more distant background higher 
above the Galactic Disc. Mapping of the nearby ISM in- 
deed suggests that the path towards ui Cen first traverses 
the warm/hot gas of the Local Bubble until ~ 100 pc, then 
passes through a neutral medium until ~ 200 pc before 
leaving the Disc (Lallement et al. 2003). Although wCen 
is not seen through the bright Ha emission of the Galac- 
tic Disc, it is seen through a layer of fainter, more diffuse 
Ha emission (Gaensler et al. 2008) possibly associated with 
an extra-planar warm ionized medium. Given that Na I and 



Can both trace the warm neutral medium one might have 
expected that the absorption bears a prominent imprint of 
the slab of neutral cloud material. 

This does not appear to be the case. The column 
density ratio Nc a iiK/NNaiD2 ~ 1.7-3.4. For warm, inter- 
cloud material, Crawford (1992b) found that Nca/NNa ~ 
0.036 NcaiiK/NN a iD2- In our case this would suggest that 
Nca/NNa ~ 0.06-0.12, i.e. a depletion of calcium by a factor 
~ 10. The ratios observed here suggest clouds of T c 7000 K 
and (probably within an order of magnitude) n c ~ 1 cm -1 , 
also on the basis of calculations by Bertin et al. (1993) , which 
is typical of the warm ionized medium. 

Wood & Bates (1993, 1994) and Bates et al. (1992) have 
probed the ISM in front of, and surrounding, u> Cen with a 
few dozen stellar probes at distances from ~ 100 pc to > 0.5 
kpc. They identified four main kinematic components which 
are also spatially distinct: (i) absorption around zero km s _1 
arising from hot gas in the Local Bubble; (ii) absorption 
around —10 km s _1 arising from cooler material between the 
Local Bubble and the expanding Scorpius-Centaurus shell; 
(hi) strong absorption around —15 km s _1 , associated with 
IR emission seen in IRAS maps, in a vertical extension of 
the Carina-Sagittarius spiral arm at ~ 1-2 kpc distance from 
the Sun, and (iv) weaker and patchier absorption at more 
negative velocities, typically —40 km s~\ arising from extra- 
planar gas beyond the Carina-Sagittarius arm, at J> 300 pc 
above the Disc. If most of the absorption seen in the Ca II K 
and Nal D2 lines is around 1-2 kpc distance then we probe 
scales of ~ 1-20 pc. 

The mean velocity of the Call K line is (wLSR.CaiiK) = 
—27 km s" 1 . This suggests that the Call K line originates 
predominantly in the extra-planar gas at a few kpc distance, 
either because it is warmer there or because dust depletion 
diminishes the Call K line closer to us. The mean velocity 
of the Nal D2 line is («LSR,Nai D2) = —14 km s _1 . The Nal 
D2 line thus appears to arise predominantly in the puffed-up 
spiral arm at 1-2 kpc distance. This picture is confirmed by 
the fact that the more negative velocity seen in the Call K 
line compared to the Nal D2 line occurs where the Call K 
line is broader, the WcaiiK/WNaiD2 ratio is higher and the 
reddening is lower. 

The kinematic lag observed in the extra-planar gas may 
not just be due to the projection of Galactic rotation. Gas 
infalling from the Halo will lag behind the Galactic rotation 
(Albert et al. 1993; Shull et al. 2009). Crawford et al. (2002) 
performed a detailed study in the Call K and Nal D2 lines, 
but at much higher spectral resolution, of similar warm gas 
in the lower Halo seen through the Local Interstellar Chim- 
ney. They found evidence of cooling gas precipitating onto 
the Disc. Alternatively, drag experienced in the Disc-Halo 
interaction zone where the rotating gaseous Disc meets the 
hot Halo gas might slow down rotating extra-planar gas. 
Such interaction might also give rise to instabilities in the 
Disc, resulting in a spiral density wave, much akin to the 
standing sound waves generated when rubbing a metal pipe. 

To investigate the possible presence of a velocity lag, 
we computed the line profile analytically, from: 

,- t (»lse) 



I(v) 



(2) 



where the optical depth is described as a function of distance 
to the Galactic Centre, r, and the Galactic plane, z, in a 
double exponential form: 
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Figure 15. Analytical line profiles for different scaleheights, and 
for the case of drag exerted by the hot Halo gas. They all peak 
sharply, either at ^lsr = —53 or at —42 km s — 1 . The profiles are 
also convolved with an instrumental profile appropriate for our 
spectra (dashed), and then peak around «lsr ~ ~~ 30 km s _1 , to 
be compared with the observed Call K profile (coarse histogram). 
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where h r and h z are the Disc scalelength and scaleheight, 
respectively, and d is the distance from the Sun. 

The analytical expression for d(«LSR) assuming a flat ro- 
tation curve, v(r) — vq, and some of the derived quantities, 
are discussed in Appendix A. There, we also describe the 
more ungainly but still analytical result for the case where 
we introduce a velocity lag as a function of z ( "drag" ) : 



v(z) =v @ (l - -^-sin(fc)) , 



(4) 



where X z sets the rate of linear decline in the velocity. We 
assumed a distance to the Galactic Centre of r© = 8 kpc, 
v @ = 250 km s" 1 (c.f. Shattow & Toeb 2009), and h r = 3 
kpc. The results are visualized in Fig. 15, compared to the 
observed Call K line. 

The model with h z = 1.8 kpc (c.f. Gaensler et al. 
(2008), for the warm ionized medium) peaks at too nega- 
tive a velocity. The peak shifts to a more positive velocity 
if the scaleheight is lowered, or if some drag is introduced. 
The negative-velocity wing of the line profile shifts more in 
the latter case. Our average spectrum of the Call K line is 
marginally more consistent with a lower scaleheight without 
the need to invoke drag. On the other hand, the most neg- 
ative component in the higher resolution spectra presented 
in Bates et al. (1992) peaks at ulsr = —42 km s _ , which is 
exactly where a model peaks that has a moderate amount of 
drag (X z ~ 30 kpc). In that model, this velocity extremum 
is reached at d = 4.8 kpc, i.e. just in front of u> Cen (assumed 
to be at d = 5 kpc) . 



4.3 Extra-planar small-scale ISM structure 

As concluded in the previous section, our view towards 
u) Cen seen through the eyes of the Ca II and Na I lines and 
A5780 and 5797 DIBs mostly shows the spiral arm-Halo in- 
terface and the extra-planar gas beyond it. Although we see 
neutral atomic and molecular g well as dust deple- 

tion, that portion of the Disc-Halo interaction region has 
a prominent component of warm ionized gas, with a scale- 
height of ~ 1.8 kpc and a filling factor of ~ 30% (Gaensler 
et al. 2008) — note that ui Cen is 1.3 kpc above the Galac- 
tic plane. The hot corona permeating the Halo, making up 
most of the remaining 70% in volume, is not probed here. 

Even though the ISM is a dynamic environment in 
which pressure equilibrium is unlikely to be strictly valid, 
it is largely the case that the volume filling factor of an ISM 
phase matches its temperature, following the equation of 
state P oc pT/fi (with jj, the mean molecular weight). Cold 
media are found in the form of compact structures, with a 
high density (exacerbated by the lower temperature). But 
as the cross-section of such structures decreases, it is harder 
to sample the colder medium by means of discrete sightlines 
(as opposed to full imaging). Furthermore, superposition of 
separate cold structures is likely if the column length is sub- 
stantial. Variations in column density across the sky should 
be reflected in variations along the sightline. 

We have performed a simplistic analytical modeling at- 
tempt to assess the small-scale structure of the extra-planar 
neutral and low-ionized medium. This assumes a size distri- 
bution of spherical clouds of radii r, between ro and Rq: 



with the specific density postulated by 



p( r ) = Po — 
\ r 



(5) 



(6) 



The normalisation of the size distribution is obtained from 



dV(r), 



(7) 



where / is the volume filling factor for a rectangular vol- 
ume of space, V. In Appendix B we describe the derivation 
from this of the expected relative fluctuations in column 
density, n. For reasonable assumptions of the parameters 
upon which this depends (see Appendix B), one would ex- 
pect (<rn/n) ~ 0.1-0.2, and this is entirely consistent with 
the fluctuations we see in our maps. This once more con- 
firms that we detect pc-scale structure in the neutral and 
low-ionized extra-planar gas. 

We can not ascertain whether the extra-planar small- 
scale structure originates in the Galactic Disc and simply 
persists as gas migrates to higher elevations, or whether it 
is generated in situ, as a result of thermal instabilities in 
cooling gas. The latter may be associated with gas raining 
down from outflows generating within the "Galactic Foun- 
tain" paradigm (Shapiro & Field 1976; Spitzer 1990) — c.f. 
Crawford et al. (2002) for recent evidence of this scenario. 
Alternatively, it may drizzle down from much further afield 
(e.g., Stanimirovic et al. 2006; Bland-Hawthorn et al. 2007). 
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5 CONCLUSIONS 

Medium-resolution spectra of 452 blue-Horizontal Branch 
stars in the globular cluster u Centauri were obtained, with 
the aim to determine the spatial and thermodynamical 
structure of the intervening diffuse interstellar medium. The 
study presented here is based on the Ca II K and Na I D2 lines 
and the diffuse interstellar bands at 5780 and 5797 A, which 
probe mainly the warm neutral and weakly-ionized medium. 
Maps were created spanning 0.7° , showing fluctuations in all 
these tracers, up to a factor two in the Ca and Na lines, on 
all scales down to a few arcminutes. These fluctuations were 
reproduced in a statistical sense adopting a simple analyti- 
cal model for a population of clouds with power-law distri- 
butions over size and density. Correlations with reddening 
confirmed that calcium depletion onto dust grains occurs in 
neutral clouds. The location of the Ca K and Na D dom- 
inant line- forming regions was deduced, on the basis also 
of kinematical evidence (line centroids and line widths) in 
agreement with literature high-resolution studies of a small 
number of sightlines. 

The most pertinent conclusions of our investigation can 
be summarised as follows: 

• In the direction towards w Centauri, (I, b) — (309°, 15°), 
the Ca K and Na D lines and the 5780 and 5797 A DIBs trace 
predominantly the neutral and weakly-ionized medium in 
the extra-planar Disc-Halo interface. The neutral medium 
is relatively more prominent in the front of this extended 
region, ~ 400 pc above the nearby (1-2 kpc) Carina- 
Sagittarius spiral arm, whilst the ionized medium dominates 
behind it, ~ 0.5-1 kpc above the Galactic plane. 

• The extra-planar gas is dominated by cr-type clouds, in 
which the macromolecules held responsible for the DIBs are 
electrically charged. The 5780 A DIB is identified with such 
carrier; it also appears stronger in the envelopes of neutral 
clouds, i.e. in the contact interfaces with warm gas, whilst 
the 5797 A DIB is stronger within neutral clouds. This is 
consistent with the 5797 A DIB being associated, either with 
a neutral carrier molecule or with a singly-charged cation 
at a lower photo-ionization threshold of the parent neutral 
molecule. 

• Small-scale structure is detected in both the neutral 
and weakly-ionized medium, down to scales of ~ 1 pc. This 
implies that small-scale structure persists — or is generated 
— well above the Galactic plane, several hundreds of pc (or 
more) away from the nearest sites of supernovae. 
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APPENDIX A: COMPUTATION OF LINE 
PROFILES WITH AND WITHOUT DRAG 

To compute a line profile (Eq. (2)), we need to assess the 
optical depth in a velocity interval. This is related to the 
column density, and thus obtained from 



t = r n(r, z) 



dd(WLSR) 



dl>LSR 



(Al) 



In our application, the precise value of to is not important 
so we scale it to match the observed line profile. The density, 
n(r, z), as a function of distance to the Galactic Centre, r, 
and Galactic plane, z, is described in our model by a double 
exponential (c.f. Eq. (3)). The main mathematical exercise 
is to prescribe r, z, and distance d as a function of wlsr- 
Trigonometry relates r and d: 



r 2 = r-Q + d 2 cos 2 (&) — 2r Q dcos(7) cos(6), 
and z and d: 
z — dsin(fc), 



(A2) 



(A3) 



where I and b are the Galactic coordinates of the sightline. 
Further simple trigonometry yields a prescription for the 
radial velocity difference as a function of r (and thus of d): 



ulsr = ( W ~~T ~~ w ©) sin(Z) cos(6), 



(A4) 



where v and vq are the Galactic rotation velocities of the 
gas cloud and the Sun, respectively. For a flat rotation curve, 
v (r) — vq, the distance at which the velocity difference 
reaches an extremum is given by 



re 



cos(Z) 
cos(b) 



(A5) 



For the sightline towards w Cen, (l,b) = (309°, 15°), and 
adopting = 8 kpc and vq — 250 km s~ , one obtains 
d x = 5.2 kpc and a corresponding «lsr,x = —54 km s _1 . 

To obtain a prescription for d(«LSR), we need to invert 
Eq. (A4) (incorporating Eq. (A2)). For the case of a fiat 
rotation curve this is straightforward: 



d(vLSR) = (cos(Z) ± (a 2 - sin 2 (Z)) 



cos(6) ' 

where we introduce for the sake of brevity: 



ulsr 



vq sin(Z) cos(&) 



+ 1 



(A6) 



(A7) 



The plus/minus sign in Eq. (A6) arises from the fact that 
beyond d x the velocity evolves back towards zero, resulting 



in two solutions for the distance at any given velocity. In our 
situation, this is avoided by us only tracing gas up to the 
distance of u> Cen, which is just before the turning point d x . 
Hence, in our case the minus sign applies to Eq. (A6). The 
derivative is also readily derived, irrespectively of the sign: 



(A8) 



The situation becomes more laborious if we replace the 
flat rotation curve by the modification following Eq. (4), to 
allow for some kind of drag between the Disc rotation and 
the Halo. However, the choice of this form still allows an 
analytical solution to be obtained, viz. for the distance: 

ro 



dd(«LSR) 




o 3 r 


d^LSR 




vq sin(Z) cos 2 (b)(a 2 - sin 2 (Z)) ' 5 



d(fLSR = (p3 ± (pi -P1P2) ' 5 ) 



and for the derivative: 



P2 cos(6) ' 



dd(«LSR.) 




d«LSR 





vq sm(Z) cos (0) p2 



(A9) 



(A10) 



where we have made use of the following abbreviations: 

(All) 



k = --2- tan(b); 

A z 



Pi = 1 - a ; 
P2 = 1 
pz = cos(Z) 



, 2 2 
k a ; 



ka 



. A: 



qi = -2— (pz ± (pi - pip?) 5 ) ; 

P2 

<?2 = 2fc ± (2kpz - k 2 p! - P2)(p 2 - piP2)~°' 5 - 
Now, the turning point will have shifted, to 
/ j cos(Z) — k \ r Q 



1 — k cos(Z) J cos(b) 



(A12) 
(A13) 
(A14) 

(A15) 

(A16) 

(A17) 



For moderate drag, with X z = 30 kpc, the turning point is 
located at d x = 4.8 kpc, and «lsr,x = —43 km s _1 . Hence 
we must add the contributions from both solutions: the one 
up to d — d x , with the minus sign in Eqs. (A9), (A15) and 
(A16), and the one from d x until 5 kpc, with the plus sign. 



APPENDIX B: DERIVATION OF EXPECTED 
COLUMN DENSITY FLUCTUATIONS 

The volume contributed by clouds (following Eqs. (5)-(7)) of 
a size between r and r+dr is dV(r) = jnr 3 dN(r). Assuming 
the volume V has a width R and depth L, we obtain 



A, 



where 



(3 - (R 3 - ( - 

(HRo/ro))- 1 



for £ 3, 
for £ = 3. 



(Bl) 



(B2) 



The expectation value for the column density, (n) , that 
we expect to measure along a given sightline is obtained by 
integrating the column density dn resulting from clouds of 
size r, which is given by 



dH = 4rp(r)(r/R) 2 dJV(r). 



(B3) 
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Here, a factor n(r/R) 2 was introduced to reflect the cloud's 
relative cross section, and a factor 4r/ir to reflect the mean 
path through a spherical cloud with radius r. Integrating 
over all cloud sizes, we thus obtain: 



<II) = 4 Po N 



_o 

R? 



flu 



„2-«-x 



dr. 



(B4) 



Fluctuations in the number of clouds of a size between 
r and r + dr along a given sightline, 

tIn =^(^) 2 diV(r), (B5) 

cause fluctuations in the column densities amongst differ- 
ent sightlines, (on). The expectation value for the variance 
in column density is obtained by adding in quadrature the 
contributions to the variance by clouds of different sizes, viz. 
((Ar/7v)p(r)a dN ) 2 : 



I 2 \ 16 2 Ar 



,.e+2x 

_0 

R 2 



„3-f-2x 



dr. 



(B6) 



where S is the separation between sightlines. The reason 
for not integrating up to r > Ro is that sightlines behave 
covariantly on a scale that is smaller than the cloud size, 
hence the contribution of clouds to the variance diminishes 
if they are larger than the typical sightline separation. We 
thus obtain for the expected relative fluctuations: 



<<Wn> = 

where 



\ULA 



1/2 



"1/2 



B, 



(B7) 



B = < 



(ln^/rp)) 1 / 2 
ln(fl /r ) 



for e j= A n # 0, 



for e = A T) = 0, 



1/2 



(B8) 



, ,„ , — r for e 7^ A n = 0, 



(ln(<5/n,)) 1/2 - 



for e = A r\ 0. 



where e = 4 — £ — 2x and 77 = 3 — £ — X- 

As an example, for (, — 2 and x = 1 ( c -f- Kim et al. 
2007), i.e. e = 77 = 0, and reasonable values for Ro ~ 10 
pc, r ~ 1 AU, 5 ~ 6 pc (10' at 2 kpcfl L ~ 2 kpc, and 
/ ~ 0.3, we obtain (an /II) ~ 0.084. It largely depends on 
(/ L/Ro) -1 ^ 2 (as long as Ro S> ro), so for a larger cut-off in 
cloud size and/or a shorter characteristic column the value 
for (crr/II) increases. As Ro < R, with R ~ 25 pc (0.7° 
at 2 kpc), and L is unlikely to be much less than a kpc, 
(<7n/n) could possibly reach ~ 0.2. It also increases if the 
filling factor is smaller, so if the contribution of colder neu- 
tral medium is relatively more important than the warm 
ionized medium then / < 0.3. This is likely to be the case, 
as we saw in the previous sections that we probe a mixture 
of neutral and low-ionized material (c.f. Welsh, Wheatley & 
Lallement 2009, who probe similar conditions). 



For n sightlines distributed randomly over an area of R 2 , the 
average separation is (5) ~ R/y/n.. 



